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Abstract 

We show that there exist asymptoticaUy flat almost-smooth initial 
data for Einstein-perfect fluid's equation that represent an isolated 
liquid-type body. By liquid-type body we mean that the fluid energy 
density has compact support and takes a strictly positive constant 
value at its boundary. By almost-smooth we mean that all initial 
data fields are smooth everywhere on the initial hypersurface except 
at the body boundary, where tangential derivatives of any order are 
continuous at that boundary. 

PACS: 04.20.Ex, 04.40.Nr, 02.30.Jr 



1 Introduction 

There is still missing a description by an initial value formulation in general 
relativity of a self-gravitating ideal body in a general situation. By an ideal 
body we mean a perfect fluid where the thermodynamical variables and the 
fluid velocity have spatially compact support. Examples are, within some 
approximation, a star, a neutron star, or a fluid planet. By a general situation 
we mean first, a body without symmetry, because spherically symmetric 



bodies are already described, whether static or in radial motion |T7 
and second, a situation including nearly static objects. 
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It is remarkable the lack of this type of description. Stars are common 
objects in the universe, and a perfect fluid is the simplest matter model 
for them. General relativity is the currently accepted theory of gravitation 
to describe stars, as well as planets, white dwarfs, and neutron stars. An 
initial value formulation will be a useful tool to predict the time evolution of 
such objects as predicted by Einstein's equation, without any approximation 
besides the choice of the matter model. 

The main difficulty is to find a solution in a neighborhood of the time-like 
hypersurface corresponding to the fiuid-vacuum interface, where Einstein- 
Euler's equation transforms into vacuum Einstein's equation. It is known 
how to describe regions not including this interface, since an initial value 
formulation for vacuum Einstein's equation was first given in and for 
Einstein-Euler's equation with everywhere non-vanishing energy density, in 
0. The problem at the interface is inherent in the fiuid equations and it is 
also present in a Newtonian description. A summary of known results on free 
boundary problems is given in Sec. 2.5 in [^], in the context of Newton's 
theory as well as of general relativity. 

A first step to set up an initial value formulation for Einstein-Euler's 
equation in a neighborhood of the fiuid-vacuum interface requires to find, 
from the complete system of equations, a symmetric hyperbolic system that 
remains symmetric hyperbolic even at that interface. A first system of this 
type was found in ||25| for a certain class of fluid state functions. However, 
spherically symmetric static situations can not be described by these state 
functions. The reason is that every solution found in that reference satisfies 
that the fiuid particles at the body boundary follows a time-like geodesic. 
But this is not the case for a spherically symmetric static star, as the fol- 
lowing argument, already given in pSf, shows. Consider a static spherically 
symmetric stellar model. The fiuid 4-velocity must be proportional to the 
timelike, hypersurface orthogonal Killing vector. The proportionality factor 
achieves that the 4-velocity be a unitary vector field. In the vacuum region 
the space-time must coincide with Schwarzschild's. At the star boundary 
the timelike Killing vector field must coincide with the timelike Killing vec- 
tor in Schwarzschild's space-time. However, the 4-velocity obtained with 
this Killing vector field is not geodesic. Therefore, spherically symmetric 
stellar models are not included among the solutions given in |2^ and so, one 
does not expect that nearly static stellar models can be described with these 
solutions. 

A second system of the type mentioned above was found, for a general 
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class of state functions, in IT^. However, the initial value formulation that 



describes nearly static objects is still missing. Fix some smooth state function 
such that both the energy density and the sound velocity vanish at pressure 
zero, and then consider the system given in [|14| for that fluid. Assume that 
there exists a smooth solution that describes such a fluid body. Then, one 
can check that the fluid particles at the boundary of the body follow a time- 
like geodesies. In other words, smooth perfect fluid solutions of system given 
in 



14 1 can not describe nearly static situations. 



Therefore, a first attempt to describe nearly static perfect fluid bodies by 
an initial value formulation would be that the fluid state function satisfies the 
following condition: Neither the fluid energy density nor the sound velocity 
must not vanish when the pressure vanishes. Because at the boundary of 
the body the pressure vanishes, we are then requiring that the border of the 
body have non-zero energy density. We call them "liquid-type" bodies, and 
"liquid-type" state functions. Therefore, it is natural to study what are the 
appropriate initial data for a liquid-type body. It turns out that the answer 
to this question was not known, it is subtle, and is the subject of this work. 

Fix once for all a simple perfect fluid, that is a perfect fluid with 1- 
dimensional manifold of fluid states; for example one with a state function 
of the form p{p), where p is the fluid pressure and p the fluid co-moving 
energy density. Assume that this state function is smooth and of liquid- 
type. By liquid-type ideal body initial data we mean a three dimensional 
initial hypersurface, its first and second fundamental forms, and the fluid 
initial 3-velocity and co-moving energy density. The first two fields must 
be asymptotically flat, the last two must have the compact support, the 
support of the fluid velocity must be included in to the support of the energy 
density; and all of them must be a solution of the constraint equations, and 
satisfy some energy condition. In addition, there exists an extra constraint 
on the initial fluid fields: The fluid co-moving energy density must be strictly 
positive and constant at the border of its support. Constant because the 
simple perfect fluid state function implies that there exists only one single 
value of the fluid co-moving energy density such that the pressure vanishes, 
and this is the value of the energy density at the border of the body. Notice 
that only the energy density as measured by a fluid co-moving observer must 
be constant. This extra condition only arises for liquid-type fluids, because 
it is trivially satisfied in the case where the fluid co-moving energy density 
vanishes at the border of the body. [See Eqs. (^)-(H)-] Also notice that 
this data must have a first fundamental form, that is, with at least one 
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continuous derivative everywhere; if not, Dirac's delta appears in the fluid 
energy density. 

Initial data of this type was not known in the literature. The only result 
on solutions of the constraint equations with discontinuous matter sources 
and first fundamental form, P|, does not guarantee that the fluid co- 
moving energy density be constant at the border of the body. Here is why. 
The solutions with discontinuous matter sources are found by the usual con- 
formal rescaling that also rescales the matter sources. Then, the initial phys- 
ical energy density [the energy density as measured by an observer at rest 
with the initial surface, function in Eq. (^] is the product of the initial 
unphysical energy density (free data) times the conformal factor at power 
minus eight. We do not know any procedure to choose the free initial data 
such that the solutions given in |^ guarantee both, that the fluid co-moving 
energy density [function p in Eq. (^)] be constant at the border of the body 
and the conformal factor be C^, simultaneously. We should mention that 
there are also found in [§ solutions of the constraints without rescaling the 
energy density, but in this case, only continuous energy densities are consid- 
ered, that is they vanish at the border of the body. 

Here, we conformally rescale all the fields except the initial physical en- 
ergy density, which is now free data, given positive and constant at the 
boundary of the body. We impose that the fluid 3-velocity vanishes at the 
body boundary. These two conditions imply that the fluid co-moving energy 
density is positive and constant at the body boundary. [See Eq. ([7p.] The 
subtle part now is to solve, with the initial physical energy density as free 
data, the equation for the conformal factor. This is a semi-linear elliptic 
equation, with the non-linear term given by a discontinuous function (the 
initial physical energy density) times the unknown (the conformal factor) at 
the power plus five. In other words, a non- decreasing function of the un- 
known, times a discontinuous given function (in contrast with [Qj where this 
function is continuous). We introduce a compact manifold, and we prove ex- 
istence of a conformal factor, based on Schauder's fixed-point Theorem. 
The proof follows the ideas given in an appendix of [^ , and in [^ . There is 
only one (technical) requirement on the initial physical energy density: Its 
L^-norm, computed with the unphysical metric, must not exceed some given 
upper bound. We show in the appendix that, although this condition ex- 
cludes possible initial data, it is mild enough to include interesting physical 
situations, such as neutron stars. 

We also give a statement on the regularity of these data. They can 
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not be smooth, because the hquid-type energy density is, by definition, a 
discontinuous function of the space variables. How regular can it be? The 
smoothest liquid-type body is almost-smooth; that is, smooth everywhere 
except at the body boundary, where tangential derivatives (appropriated 
defined) of any order are continuous. In other words, we prove the following: 
If the unphysical metric is a smooth field on the initial hypersurface, and 
the fiuid free initial data are smooth up to the body boundary with every 
derivative tangential to the boundary continuous through that boundary, 
then the same holds for all the initial data fields. A crucial requirement to 
prove this statement is that the conformal factor be on a neighborhood 
of the body boundary. (See Sec. ^]3|.) 

Summarizing, we prove that almost-smooth initial data representing liquid- 
type simple perfect fluid bodies can be obtained as solutions of the constraint 
equations, in the case where the fluid 3-velocity vanishes at the body-boundary, 
through a suitable modification of the usual conformal rescaling techniques. 

Finally, some technical remarks: (i) We rescale the initial fiuid momentum 
density but not the initial fiuid energy density, in order to solve the constraint 
equations. Therefore, we have to choose the rescaled momentum density 
small enough, in the sense given in Sec. [4.1| , in order to have physical data 
satisfying the dominant energy condition, (ii) We impose that the fiuid 3- 
velocity vanish at the body boundary. There is no physical justification 
for this assumption, it is made because is the only way we know, with the 
rescaling of the initial data field that we have chosen, to guarantee that the 
fiuid co-moving physical energy density be constant at the body boundary. 
[See Eq. (0).] We give an interpretation of this condition in Sec. ^]2[ (iii) 
Besides the conformal rescaling to solve the constraint equations we perform 
a conformal compactification in order to solve elliptic equations on some 
unphysical compact manifold. Asymptotic decay properties of fields in the 
physical initial hypersurface are translated into differentiability properties of 
these fields at a particular point in the unphysical compact manifold. These 
differentiability properties at the point at infinity are completely independent 
of the differentiability of the fields near the body boundary. 

In Sec. we introduce the main definitions we need in order to present 
the principal result. Theorems Q and 0. We also give a proof based on results 
obtained in Sees. ^ and ^ In Sees. and |3.2| we give the main existence 
proofs for the semi-linear and linear elliptic equation associated with the 
Hamiltonian and momentum constraint, respectively. In Sec. p.3| we prove 
the regularity statements. Theorems ^ and ^. In Sect O we explain why this 
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discussion on an energy condition appears, and we give a simple condition 
on the free data such that the physical initial data satisfies the dominant 
energy condition. The constraint equations are naturally written in terms 
of the initial fiuid energy density and the initial fiuid momentum density. 
Eqs. ®-(||) relate them to the fiuid pressure, co- moving energy density, 
and 3-velocity. In Sec. [4.2|, we prove that these equations are invertible. 



In Sec. ^ we comment on the initial value formulation for liquid-type ideal 
bodies. In Theorem ^ we need to assume that the fiuid initial energy density 
satisfies an inequality [Eq. (|2^) ] involving both the unphysical manifold and 
the unphysical rescaled metric. In Appendix |^, we study this inequality. In 
Sec. A.l , we show (Lemma P that a similar (but weaker) inequality holds 
for every initial data. By an explicit example, in Sect |A.2| , we prove that the 
inequality required for the existence theorems is in fact a restriction on the 
allowed initial data. This example also suggests that this restriction is mild, 
in the sense that interesting physical systems, like neutron stars, satisfy it. 



2 Definitions and main result 
2.1 Liquid-type ideal body data 

We first introduce what we mean by initial data for a liquid-type ideal body. 
Afterwards, we split the concept of almost-smooth into two pieces. Given a 
field and an open bounded set Q on some manifold, we introduce the concept 
of an f2-piecewise smooth and f2-tangentially smooth field. Finally, in the 
next subsection, we present our main result. 

Consider an initial data set for Einstein's equation with matter. That is, 
consider a 3-dimensional, smooth, connected manifold M, a positive definite 
metric, Qab, and a symmetric tensor field, p°'^, on M, together with a vector 
field, and a positive scalar function, ft, subject to the condition jaj" < fi^, 
and solution on M of 

DaP''' - = (1) 

R + {pa"? - PabP"' = 2Kfx, (2) 

where Da and R are the Levi-Civita connection and the Ricci scalar associ- 
ated with (Jab, and k, = 8n. Indices on tensors with "tilde" are raised and 
lowered with q"'^ and gab, respectively, where QacQ'^'' = 5a^- Latin letters a, 
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b, c, represent abstract indices. The fields solving (|T|)-(|^) have a meaning 
as part of a 4-dimensional space-time solution of Einstein's equation with 
matter sources. The manifold M represents a three dimensional space-like 
hypersurface such that (jab and are their first and second fundamental 
forms. This hypersurface will be a maximal slice if and only if pa"' = 0. 
The fields ft and j"" represent the normal-normal and the (negative) normal- 
parallel components to M of the stress-energy tensor. The dominant energy 
condition on the stress-energy tensor in the space-time implies jaj" < /i^ on 
M. This condition is the reason why one does not, in general, pick any qab 
and p"* and then define j" and /i by (|I])-(^. Because if one does that, the 
resulting fields form an initial data set iff the energy condition is satisfied by 
these j" and fl. It is an open question whether there exists a procedure to 
find appropriate qat and besides the conformal rescaling one. 

The initial data set is asymptotically fiat if the complement of a compact 
set in M can be mapped by a coordinate system diffeomorphically onto 
the complement of a closed ball in such that we have in these coordinates 

g,, = (l + 2m/f)5,,+0(f-2), (3) 



= oir'), (4) 

as f := bijX^x^ — > oo, where m is a constant that represents the ADM mass 
of the data. Latin letters j, /c, denote coordinates indices and take values 
1, 2, 3, while bij = diag(l, 1, 1). 

Fix as matter source a simple perfect fiuid. That is, first, introduce on 
M a non-negative scalar field p, interpreted as the fiuid co-moving energy 
density, a vector field v", interpreted as the fiuid initial 3- velocity, and fix a 
function p{p), the state function, interpreted as the fiuid pressure as function 
of the co-moving energy density. Second, introduce on M the equations 

^ (p + pK . . 

with v'^ = VaV" < 1. In the space-time solution of Einstein's equation with 
matter sources, the normal-normal and the (negative) normal-parallel com- 
ponents to M of the usual perfect fiuid stress-energy tensor are precisely 
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the left hand side of Eqs. (0)-®, respectively. If u"" denotes the unit 
fluid 4-velocity, and the unit normal to the initial hypersurface, then 
u"" = (ri" + v"') / a/1 — -y^. For this matter model the dominant energy con- 
dition is equivalent to p > p. In addition one can prove that p > p implies 
jaj"' < jj^- The sketch of the proof is the following: from Eqs. d^)-® define 
f{p,v) := VXi"//^- Notice that /(p, 0) = 0, and /(p, 1) = 1. One can prove 
that p > p implies df/dv > 0, for < v < 1; then, follows that / < 1 for 
0<v <1. 

Consider a liquid- type ideal body data set. That is, an asymptotically 
flat initial data with a simple perfect fluid whose state function is of liquid- 
type, and both the fluid 3-velocity and the co-moving energy density have 
the same compact support, n C M. By a liquid-type state function we 
mean a non-negative, non-decreasing, smooth function p{p) that vanishes 
at po > 0. As an example consider a big water drop (in order to neglect 
surface tension effects), or a fluid planet, or a neutron star. The support of 
the energy density represents the place occupied by the body. The value of 
the co-moving energy density at the border is determined by the function of 
state as the value where the pressure vanishes. (Otherwise, the acceleration 
of fluid particles lying on this border becomes inflnite.) Therefore, a liquid- 
type ideal body satisfles p\dQ = po. Eqs. (|^)-(|^) translate this condition for 
the co-moving energy density into a constraint on the flelds p and at dfl, 
where they are not longer free but they must satisfy 

[H^-lf/mon = Po- (7) 
As a summary, we state the following: 

Definition 1. A liquid-type ideal body initial data set consists of fields Qab, 
pab ^ p ^ state function p{p), such that: (i) (jab is a Rie- 

mannian metric, p""^ is a symmetric tensor, and both are asymptotically flat; 
(a) p{p) is liquid-type, and vanishes at po > 0; (Hi) suppiy"") C supp{p) = Q 
compact; (iv) p|an = po; (v) These fields are solutions of Eqs. ^-(^, ^-(^ 
on M. 

Given an open set Vl' C M"^, we denote by C^{VL') and C*'"(fi') the spaces 
of s-times continuously and Holder continuously differentiable functions, re- 
spectively, with s > integer, and < a < 1. We use the notation 
C"(fi') = C°'"(fi')- We also denote by L^iVt'), W''P{Vt'), and by Lf^,(fi'), 
^loc (^') Lebesgue and Sobolev spaces, and the local Lebesgue and local 
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Sobolev spaces, respectively, where 1 < p < oo. We follow the definitions 
given in and the generalizations for smooth manifolds, M', given in 

0. Finally, we say that a tensor field on Q' C M' belongs to one of the 
functional spaces mentioned above, if all its components, in some smooth 
atlas of M', belong to such a space. 

It is convenient to split the concept of an almost-smooth field, presented 
in the introduction, into the following two definitions. The first one is an 
fi-piecewise smooth field. Consider a smooth manifold M' , a tensor field u 
on that manifold, and a open set Q C M', with compact closure. We say 
that u is ll-piecewise smooth, if m e C°°(n) n C^{M' \ fi). Note that this 
definition involve conditions on the field both in Q and its complement. An 
example of an f2-piecewise smooth but not smooth function is any / such 
that, < / e C°°(n) and / = on M' \ H. The fiuid energy density of a 
liquid-type body is such a function. 

The second concept is an f2-tangentially smooth field. Let qab a smooth, 
positive definite, metric on M'. Assume that dQ is a smooth submanifold of 
codimension one. Let a normal vector to dfl with respect to qab- Con- 
sider a Gaussian normal foliation in a neighborhood of dQ, that is, a foliation 
orthogonal to the geodesies tangent to n"" at every point of dQ. Define •h"' out- 
side dfl to be tangent to these geodesies. Let Vqq any smooth tangent vector 
field to this foliation, i.e.; any smooth vector field such that Vg^ha = 0. We 
say that an i7-piecewise smooth field u is il-tangentially smooth if for all A; > 1 
the tangential derivatives Vq^{u) are continuous; where Vq^{u) := Vo^DaU, 
and Vq^{u) := Vg^DalVg^ for /c > 1. For example, choose the field 

to be the energy density of a liquid-type ideal body, and Q the interior of its 
support. A necessary condition for this field to be fi-tangentially smooth is 
to be constant at dQ. 

2.2 Main result 

The strategy is, first, to find fields qab, p"*, j", and /i, solution of ([l|)-(0) 
with the desired properties. Conformal rescaling techniques are used in this 
part. We also introduce a compact manifold where equations associated with 
(El)"® solved with boundary conditions chosen in such a way that, the 
de-compactification of these solutions gives asymptotically fiat initial data. 
Then we prove that under specific assumptions on the state function, Eqs. 
(H)"® '^^'^ be inverted for all p > Po and for ?)" with < v < 1. 
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Fix a 3-dimensional, orientable, connected, compact, smooth manifold. 
Fix i G M, and M := M \ {i}. The choice M = S^, and so M = de- 
scribes, for example, ordinary stars. A restriction of this type in the topology, 
however, plays no role in what follows. Let hab G C°°{M) be a Riemannian 
metric on M. Let x* be its associated Riemann normal coordinate system at 
i, and r the geodesic distance. Let hab ^ C°°{M) a symmetric tensor such 
that 

x'hij = 0. (8) 

Latin indices i,j,k denote tensor components on coordinates x\ Let qab G 
C°°{M), be a Riemannian metric on M with scalar curvature R. Let be 
an open ball of geodesic radius e centered at i. Assume that there exists 
e > 0, such that the metric qab on has the form 

qij = hij + r^hij, (9) 

in the coordinates x\ Since we have assumed (||), these coordinates are also 
normal coordinates of the metric qab- The motivation for Eq. is given in 
the remarks below Theorem 

Fix a non negative scalar field jl and a vector field on M with supp(j") C 
supp(/i) = Q, where fl is some open set with compact closure fl C M, such 
that its boundary dQ is a smooth submanifold of codimension one. Introduce 
on M the fields 6 and p"^, with pa"' = 0, solutions of 

DaP" = -k/, (10) 

= - i^C'. (11) 

where Lq{9) := q"'''DaDb9 — R9/8, and Z^a is the Levi-Civita connection 
associated to qab- Indices of "non-tilde" tensors are raised and lowered with 
and qab, respectively, where qacl'^'' = Sj'- Fix the boundary condition 

P'' = 0{r-'), (12) 
limr^ = l. (13) 

r— >0 

The main part of this work is to prove existence of solution to Eqs. (p!0|) - 
([T3|), and then to show that if the source function /i and are f2-piecewise 
and tangentially smooth, then so are the solutions and 6. Once these fields 
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9 and p""^ are known the initial data set is given by the following conformal 
rescaling 



qab = e\ab. p'"' = e-'Y\ f = e-^^f. (i4) 

Notice that we do not rescale the energy density /2, but we do rescale the 
momentum density j". In this way we achieve both that /i be a free data, 
and that the momentum constraint decouples the Hamiltonian constraint, 



respectively. One can check that if 6 and p are solutions of Eqs. ([TO|)-(pJ|) 
then the rescaled fields in Eqs. (0) satisfy Eqs. One can also check 

that the boundary conditions ([T2|)-(p!3D on 9 and p"'^ imply that the rescaled 
initial data is asymptotically flat. (See [||, |12|.) 

Let 7 be the Green function of the operator Lq given in Eq. ([TT|) , which 
is defined in Eq. (PU|)-(PT]). Our first main theorem is concerned with the 



momentum constraint (|T0]), (|T2D. 



Theorem 1. Fix M, M, Q, and qab as above. Let s"'' a symmetric trace-free 
tensor m W^'^i^M) n C~(M), g > 3. Let p""^ he given by Assume that 

(i) supp{f) C n C M. 

(ii) G L^{M) and it is Q-piecewise and Q-tangentially smooth. 
(Hi) Condition ^J^ ) is satisfied. 

Then, there exist a unique tensor p"'^ given by solution of Eq. l\TI\), 



(\T^). Moreover, p"-^ is Q-piecewise and Q-tangentially smooth and satisfies 

PabP^'h' e L\M). (15) 

The existence part of this theorem is essentially the standard York split- 
ting (cf. IQ) adapted to our setting. It is given in Theorem |, under a 
weaker hypothesis, s"''' is free data related to the arbitrary amount of gravi- 
tational radiation that can be added to the system keeping the matter sources 
fixed, p"* contains the linear and angular momentum of the initial data, it 
can also be prescribed freely unless there are conformal symmetries. In this 
case it has to satisfies condition (iii), which is the corresponding Fredholm 
condition (see the remark after Theorem]^ for a physical interpretation). The 
regularity part of the theorem is proved in Sec. [X^ . 

We have chosen the unphysical metric, qab, smooth on M \ {i}. This is 
a reasonable physical assumption. However to also impose smoothness at i 
it is too restrictive. In this case initial data for stationary space-times are 



ruled out (see |11]). The differentiability at i of the unphysical metric is 
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related with decay at infinity of the associated physical metric qab-, imposing 
smoothness at i means a restriction in the fall off which is, in particular, 
incompatible with the stationary solutions. In order to include these data, 
we have made the assumption (H). Although the functions hij and hij are 
smooth, the metric belongs to C^'"(i?e) but it does not belong to C^{B^). 
The data for stationary space-times have precisely this form (see fllj). In 
order to prove the theorem one certainly does not need smoothness of h^j 
and hij, only a finite number of derivatives. But the important point is that 
in order to prove the last part, Eq. (|15D , it is not enough to require that, for 
example, qab G C^'°'(M). Eq. ( [T5| ) is essential in order to prove our second 
theorem. We prove ([T5|) in Theorem |^. 

Proof. The metric given by Eq. (||) satisfies that qab G Vr'^'^(M), p > 3/2. 
By assumption j"" and s"^ belong to W^''^{M). Therefore Assumption (iii) 
and Theorem |^ imply that there exists p"^ G W^''^ (M), 1 < q' < 3/2, given 
by Eq. (|4l|) which solves Eqs. (10), ([121) . The hypothesis on the metric given 



in (^ and Theorem pimply Eq. ([T5|). Assumptions (i), (ii), and Theorems ^ 
and I imply that p"'^ is fi-piecewise smooth and f2-tangentially smooth. □ 

In order to write the next theorem we need to define some constants. 
Set Cp = ||PafeP"V(87^)l|L2(M), 7+ = maxfy7, 7. = mm-^ij) and j+ = 
maXfjA/jaj'^ . Let K, k, be the positive constants defined in Sec. |3]l|. They 
essentially depend on the metric qab and the manifold M. Finally, let cq > 
solution of the following equation 

where is the volume of fl with respect to the metric qab- There always 
exists a unique positive solution to Eq. ([T6|), since for eo > the right-hand 
side of (|16D is a positive, decreasing function of eo which goes to zero at 
infinity. 

Theorem 2. Assume that the hypothesis of Theorem holds. Let p°-^ he 
the tensor field given in that Theorem. Assume that R > 0. Fix a smooth 
liquid-type state function, p{p), with zero-pressure energy density, po > 0, 
compatible with condition (iii). Assume that < dp/dp < 1. Let jl be such 
that 

(i) supp{ix) = Q. fl\an = Po, and f\dn = 0. 

(ii) jl is Q-piecewise and Q-tangentially smooth. 
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(in) For e < eo, /i satisfies 



-8 K 

ej+7_ <po<fi< lon/o/.. , (17) 



Then there exist a positive solution 9 G C^'"(M) of Eqs. (|77y, ( |7^ wt/i 
sources given by jl and ep"^, where < e < eo and < a < 1. 

Moreover, the initial data computed with qah = O^qab, p"''^ = O^^^ep"''^ , 
= , and jx is of liquid-type, as stated in Definition They are VL- 

piecewise smooth, and fl-tangentially smooth, qab G C"^'"(M), p"^^ G C"(M). 
The fluid 3-velocity, v"" , vanishes at dVL. 

We use a non- typical conformal rescaling given in Eq. (0). The positive 
outcome is that, in this way, fi is essentially free data, and so, we can choose 
it constant at dVL. A negative outcome is that this fx must satisfy the bound 
([T7|). The upper bound in ([T7|) is related with the existence of solution, 
given in Theorem ^. In the Appendix we give arguments to show that this 
bound is only technical, that is, there exist solutions which do not satisfy 
it. However, the example presented there suggest that this bound will be 
satisfied for every realistic star. The lower bound in (0) is related to the 
energy condition. It is sufficient condition for the dominant energy condition 



to hold, see Sec. HTT 



A second negative outcome is that, in order to satisfy the liquid- type 
constraint (0), we impose j°'\an = 0, this implies v°'\dn = 0. In order to 
understand the implications of this condition on the motion of the fluid, 
assume that we have a simple, liquid type, fluid solution of Einstein-Euler's 
equation. That is, a 4-dimensional Lorentzian metric Qab and a unit time 
like vector field m", representing the fiuid 4-velocity, solutions of Einstein- 
Euler's equation. The boundary B of the fiuid is the 3- dimensional, time 
like, hypersurface where p = 0. Since we have a simple fiuid, this implies 
that p is constant on B, hence the vector defined by N°- = g°'^'VbP is normal 
to B, where V;, is the covariant derivative with respect to gab- By assumption 
A^" is not zero on B. Fix an arbitrary space like foliation, with normal vector 
n"; let M a member of this foliation. Define dQ = M H B, we will assume 
that both B and dfl are smooth submanifolds. The 3-velocity v"", defined by 
u"" = (n" + ■u°)/\/l — v^, will vanish at dil if and only if the following Eqs. 
hold 

Nau'^lan = 0, (18) 
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Na^'^Un = 0, (19) 

where a;" = e"'^'^'^UbVcUd is the twist of {eabcd is the volume element of 
Qab and the indexes are moved with gab)- Eq. (|T^) is a condition on the 
foliation: the slice M has to be tangent to A^"^. Eq. (|19]) is a condition on u"", 
independent of the foliation: the normal component, with respect to the fluid 
boundary, of the twist of u'^ must vanish on dfl. Eq. (^) is a consequence 
of the Frobenius's theorem (see for example |^l|) and the fact that dQ is a 



smooth submanifold and N'^ is hypersurface orthogonal. Note that cu" itself 
can be different from zero at dQ. Condition (|19]) is not time-propagated by 
m". This condition is imposed only on the initial slice, not in the subsequent 
evolution. Although it is a restriction, it is not clear if it is a strong restriction 
or not. 

Another outcome of this particular conformal rescaling is a lack of unique- 
ness of solutions to (|T])-(§|) in terms of the free data. 

We do not require that Q be connected. A non-connected domain can 
describe several compact bodies. 

This is not the most general result one can obtain with these methods. 
One can also find solutions which are not piecewise smooth, but with some 
finite differentiability in the interior of the support of p. One can even obtain 
solutions where the support of p itself has some finite differentiability. The 
obtainment of such more general data from the techniques used to get our 
result does not present a substantial difficulty but only a greater level of 
technical complication, which could obscure the main ideas necessary to find 
these type of data. 

Proof. The upper bound on fl given by Eq. (|1^) and Theorem ^ imply that 
there exists a strictly positive solution 6' = 7 -|- G C^'"(M) of Eqs. (|ll]), 
([TB|). Hypothesis (i), (ii) and Theorems |^ and ^ imply that 6 is fi-piecewise 
smooth and f2-tangentially smooth. 

Let Qab, P"'^, 3°' be as stated in Theorem |^. Those fields are also f2- 
piecewise smooth and f2-tangentially smooth, and they satisfy qab £ C^'"(M), 
pab ^ (7"(M). The lower bound in Eq. (|T^) and Lemma ^ imply that 
the dominant energy condition is satisfied, that is jaj" < p?. Assumption 
(i) implies that the liquid-type constraint (|^) is trivially satisfied. Finally, 
Theorem ^ implies that Eqs. (^-(^ are invertible. The state function p{p) 
is a smooth function of p, so Eqs. (§)-(|^) imply that the fields ?)" and p are 
f2-piecewise smooth and fi-tangentially smooth. Eq. (|^) and Assumption (i) 
imply that v'^\dQ. = ^- D 
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3 Existence and regularity 



3.1 Hamiltonian constraint 



Consider Eqs. ([TT|), (p^. To obtain a solution 6 we first transform tfiis 
problem on M with a singular boundary condition at i G M, into a regular 
problem on M for another function. The metric Qab has strictly positive 
scalar of curvature R and the assumption given in Eq. (|^) implies that 



Qab G ly^'^(M), p > 3/2. Therefore, Lemma 3.2 and Corollary 3.3 in |jT2 



imply that there exist a unique, positive solution 7 G C^'°(M) of the equation 

L,(7) = -47r5„ (20) 

where 5j is Dirac's delta distribution with support at i. It is also true that 

7-1 e C"(M) and 

limr7 = 1. (21) 

We introduce the function = 6 — 'y. Then, Eq. ([TT|) for 6 on M becomes 
the following equation for on M, 



Before stating the Theorem concerning existence of solutions to Eq. (E^p, 
we need some notation. Given any function g G W'^''^{M) and the operator 
Lq, we introduce k to be the constant such that \g\cO{M) < ^||-^<?(fi')||L2(M)- 
This constant can be written as = CgCi, where the Sobolev coefficient Cg is 
the constant such that \g\co{M) < Cg ||fi'||H/2,2(jv/), while cl is the constant of the 
elliptic estimate ||5'||iy2,2(7(^) < cl \\Lg{g)\\i2(^M)- (See [|1])|1]-) We introduce, as 
well, the constants Cp := WPabP''^ / {^i'^)\\l'2{m), and 7+ := supfy7. Therefore, 

PabP'^'h' e L2(M) (23) 
is equivalent to the condition Cp < 00. 

Theorem 3. (Existence) Let M and M be as in Sec. Let qat be a 

Riemannian metric on M , such that qat G l^^'^(M), p > 3/2, and R > 0. 
Assume thatp"-^ satisfies that Cp < 00. Let jibe a positive function of compact 
support in Q G M, such that 
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where K = 4P/ {5^K,k) . Then there exists a non-negative solution d G W^^'^(M ) 
of equation ^^). The solution is strictly positive unless both p""^ and jl are 
zero. Moreover, it satisfies d < (7+ + 5kCp)/A. 

Remark: The proof is based on Schauder's fixed-point Theorem (see for 
example Let B G X be a nonempty, closed, convex set in a Banach 

space X, and F : B B be a continuous mapping. If F{B) is precompact, 
then F has a fix point. The construction of the functional F is similar to the 
one made in ||12| for Theorem 3.4. The only difference is the choice of the set 
B, and the main work is to prove that for this choice we have F{B) C B. 

Proof. Consider X = C^{M), which is a Banach space under the supremum 
norm. Given a constant c > 0, define B^ = {u ^ X : < m < c}. One 
can check that B^. is a convex and closed. Define a non-linear operator 
F : Be ^ X , hj setting 

F ■■= o / 

where the f : Be ^ L'^{M) is the continuous map given by 

Under the assumptions qat G ly^'^(M), p > 3/2, and i? > it has been 
proved in ||12| that the non-linear map F is continuous and F{Bc) is pre- 



compact. The only difference between the map F and the analogous map T 
defined in |jT^ is the second term in the right hand side of (^51). This term is 
continuous. Note that 7 is singular at i but we assume that jl has support 
in fl and the point i is not included in Q. 

We only have to choose the constant "c" such that F{Bc) C Be. The rest 
of the proof shows how to find "c". In what follows we will use Lemma 3.1 
of |T^ many times. 



Introduce the functions (fc G L'^{M) and 0c ^ W^''^{M) as follows: 



Then, for all u & Be we have that f{u) — ^Pc > 0. This is equivalent to 
Lq{F{u) — (pe) > 0, and then F{u) < (pe- We now choose the best constant 
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c" such that (j)c < c. This is done as follows. Given the bound 



< k\\Lgi(i)c)\\mM), 

< k 



Cp+^(7+ + c)'||/2||L2(f,) (26) 



we impose that the right hand side of (^6[) be less or equal to c. We then 
obtain 

4 c-kCp 

This inequality has to be valid for some c, in particular for its maximum 
value given by cq = (7+ + 5A;Cp)/4. Eq. ( ]27|) evaluated at cq gives Eq. (p4D. 
Therefore, choosing B = B^o, condition (^) implies F{B) C B. Finally, 
Schauder's fixed-point Theorem implies that F has a fix point in B. This fix 
point is the solution {}. □ 

We now show that, under slightly stronger assumptions on the source 
functions /t and p"''', the function 6 belongs to C^'°(M). (This differentiability 
is important for Theorem ^.) 

Theorem 4. [C^'"(M)-regularity] Assume the hypothesis on Theorem 
hold, and let 6 = ^ + 1!}, with G ly^'^(M) solution of In addition, 

assume that jl G L'^iP), with q > 3, and that PabP°'''/l'^ G ^od^)- 
Then, 9 G Wl^{M) C ^^'"(M) zs a solution of Eqs. gl 



Proof. From the hypothesis on fi and p""^, we have /("f?) G L^j^^(M). Elliptic 
regularity implies 'd G iy^'''(M). (See fl^.) Sobolev embedding and g > 3 
imply ^ G Ci'°(M). Therefore, 7 G Ci'"(M) imply that 6 G Ci'°(M). Eq. 
(|2l| ) implies that 9 satisfies the boundary condition (plSf ). □ 



3.2 Momentum constraint 



Consider Eqs. (|Ty), (|T2|)- The main idea is, as in Sec. |3.1| , to transform these 



equations on M with a singular boundary condition into an equation on M 
for a regular variable. Solutions of this regular equation can be found by 
the transverse, traceless decomposition of symmetric tensors. See |^ for a 



transverse decomposition, and for a transverse, traceless decomposition. 



See also 0, and references therein. 
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All this procedure is performed, however, not in Eq. (|TU|) itself, but in 
a properly conformal rescaled version of that equation. The new rescaled 
metric is chosen such that its Ricci tensor vanishes at i. (The restriction 
that the unphysical metric, Qab, have strictly positive Ricci scalar on M is 
not needed in this subsection.) The positive outcome of this new rescaling is 
that it is not hard to prove that solutions p"* with non-vanishing total linear 
momentum are included. 

The plan of this subsection is, first, to introduce some notation; second, 
to set up the procedure to prove existence of solutions p""^ in a weak sense 
(Theorem and third, to prove that, under a slightly stronger assumption 
on the matter source, the solution satisfies Eq. (|23|) . (Theorem ^j.) 

We start with the new conformal rescaling. Let M, M, qab, x^, r, and 



-Be as in Sec. 2.2. Let x be a cut function, that is a smooth function with 



support in i?2e and such that x = 1 in 5^. Fix on M the metric Qab given by 

Qab = 9ab, (28) 
where the conformal factor cuq has the form 

^o = e>^^°, fo = ^x^x''L,k{t), (29) 
and we have evaluated at i the tensor field 

Lab '■= Rab — -^Rlab, (30) 

with Rab the Ricci tensor of qab- Therefore, qab = qab on M \ i?2e, and they 
differ only on B2e- One can check that Rabc^if) = 0, that is the Riemann 
tensor of qab evaluated at i vanishes. (An explicit computation shows Rabif) = 
0. Since qab is a 3-dimensional metric, Rabc'^ii) = 0.) This property implies 
that in its associated Riemann normal coordinate system at i, x-', the metric 
qab has the form 

g,, =(5,,+0(f3), f,^ = 0(f2), (31) 

where f is the geodesic distance from i measured by qab- This is the reason 
for doing the new conformal rescaling. 

We complete the rescaling introducing the fields p""^ and j° as 
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Therefore, Eqs. (PH]), (|T^ transform into 

DaP''' = -f^j\ p^^=0(f-^), (32) 

where Da is the metric connection associated to qab- Latin indices on "hatted" 
quantities represent components in the coordinate system 

We now start the procedure to transform Eq. (|32D on M with a singular 
boundary condition, into an equation on M for a regular variable. The 
singular behavior at i of a solution 9 of Eqs. ([Tl|) , (|T3|) was captured by the 
Green function 7. In the case of Eq. ([3^), the role analogous to 7 is played 
by a tensor p"-'^. The construction of this tensor field, that follows, is detailed 
in []T^, Sees. 4.1-4.2, but we briefly sketch it here. 

Consider the manifold (M, qab)- Let i?2e C U a. ball of g^b-geodesic radius 
2e centered at i, and x the associated cut function, that is, a smooth function 
that vanishes on M \ i?2e and x = 1 in B^. Let fia = DaV. Introduce on 



\ {i} the tensor field^ 

'^(1) = [SQ^'^n") - (r^ - h-h')n,Q'] , (33) 

'^(2) = ^('^"'-3^"^')' (34) 

0g) = ^n(»e'')^'^J,n,, (35) 

0^^^^ = — ^ [2P("n'') + (5'*'' - 5n"n'')neP'=] , (36) 



where A is constant, and P*^, J*^, and are constants in the coordinate 
system £*. Here = UaS""^, and in Riemann normal coordinates, = /r. 
These tensors are transverse and traceless with respect to the fiat metric. 
Let p^^^ := x(0^i;) - q^''qcd(l)lt)/^)- Finally, introduce p""^ as follows 

4 

k=l 

By construction, the tensor p"^ depends on 10 parameters, is smooth on M, 
vanishes on M\B2i, is symmetric and gab-traceless, and satisfies p*-' = 0{r~'^), 

^In the particular case where qab ~ Sat, that is the flat metric, these tensor fields 
arise as appropriate derivatives of the tensor 7^'' :— {76 J' + fiafi^) / {8f) which satisfies 
(Lgjb)"- = {-ATTSi)Sb°-. The operator Ls is defined after Eq. (Q. 
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as f — s> 0. It also satisfies 

Z)„p«'' = 0{f-^) as f ^ 0. (38) 

The last equation is obtained as follows: write DaP"'^ explicitly, and then note 
that first, the tensor fields (pl^^^ are divergence and trace free with respect to 
the flat metric, and second, that in the coordinates the metric connection 
coefficients satisfy (|3T|). 

We finally recall some needed properties of conformal Killing vector fields. 
Consequently, this paragraph is applicable to both qab and (jab- We point out 
the differentiability of the various fields, for later purposes. Fix a manifold 
{M,qab), with Qab e ^^'^(M), with p > 3/2. A conformal Killing vector 
field, is defined by (CgO'"^ := 2[D('^^^) - q'^^D^^yS] = 0, where Cg is the 
conformal Killing operator associated to the metric Qab- There are at most 
ten conformal Killing vector fields for a 3-dimensional metric. Given a vector 
field u"' G {M), with p' > 1, we say that it is orthogonal to if 

e,^" dV = 0, (39) 

M 

where the volume element is computed with the unphysical metric Qab- Notice 
that the differentiability assumption on the metric implies that ^" G C^'"(M). 
This, in turn, with the Holder inequality, implies that the integral above is 
well defined. We also introduce the conformal Killing data at i, that is, 

ka = ]:D,Dbei^, S'' = e'^'^DbU^, = Cit), a = l-DaCi^- (40) 
D 3 

Since M is connected, the integrability conditions for conformal Killing fields 
(cf. [0) entail that these ten "conformal Killing data at i" determine the 
field uniquely on M. 

We have the following existence theorem, which is a generalization of 
Theorem 16 proved in |1T2|. 



Theorem 5. (Existence) Let M, and M be as in Sec. 2^. Assume qab G 
W^'P{M), p > 3/2. Let p"^ he defined by and and qab as in Let 

gab g |yi:P (^M^ 5e a symmetric traceless tensor, and G (M), with 
p' > 1. 

(i) If the metric qab admits no conformal Killing vectors on M, then 
there exists a unique vector field G W'^''^{M), with q = p' if p' <?)j1 and 
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l<g<3/2 if p' >?)/2 such that the tensor field 

P<^b = ,^10 [p-b + S"'' + (CgW)"''] (41) 



satisfies Eqs. (17^, (III) 



(a) If the metric qab admits a conformal Killing vector i^" on M , corre- 
sponding to the conformal Killing data given in Eq. ^^), then a vector field 
as specified above exists if and only if the following condition holds, 

P''ka + JaS'' + Aa + {P' W'ii) + Q'')qa = n [ jaC dV, (42) 
where the constants P", Ja, A, and Q° characterize the tensor as in 



Proof. Because of ( |38| ) we can consider Daij)""^ + s""^) as a function in L'^{M), 
1 < g < 3/2. The equation Dap""^ = —nj^ is equivalent to 

Dair' + + {'C^wT'] = -KU^'^f, (43) 
which can be written like 

{Lgwf = -Ku^^^^j' - DaiP"' + s"'), (44) 

where (L^w)"' := DfjlCgw)""^ is an elliptic operator. Its kernel consists of all 
conformal Killing vectors, of qab, and so, of qab- Following |12[, one can 
prove that the right hand side of (^^ is orthogonal [in the sense given in 
(PPj)] to every conformal Killing vector field, if and only if Eq. (^5]) holds. 
Therefore, the assumptions on the metric, qab, and the Fredholm alternative 
for this operator imply there exists a unique solution w"' G W^^'^(M). (For 
smooth metric this is a standard result, for metric in the Sobolev space 
W^'P{M) see §.) □ 

The quantities and J" in tensor p"-'^ represent the total linear and an- 
gular momentum of the data. These quantities can be prescribed freely in 
case (i), so they are not related with the matter sources j". The interpre- 
tation is that gravitational waves can carry an arbitrary amount of linear 
and angular momentum. In the case that the unphysical metric has confor- 
mal symmetries these quantities are restricted by condition (^2]). In order to 
understand this condition, consider the case where only one Killing vector 
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exists, and it is a rotation. That is, only S"" is different from zero. We 
can always choose S"" to be a unit vector. (This vector is parallel to the 
axis of the rotational symmetry.) Construct the following initial data: first, 
choose any J"- pointing in the same direction as 5", and second, choose the 
other part of the free data preserving the symmetry. Then, all the fields in 
the initial data set have this symmetry, and therefore the whole space-time 
obtained from this initial data set will also have a Killing vector suitably 
extended outside the initial hypersurface. Condition (^21) reduces to 



J = n jaC dV, (45) 



M 



where J = a/ JaJ"^- Eq- ( ^5)) is just the standard Komar integral. (See 
for example ||3l]].) This is consistent with the interpretation that axially 
symmetric gravitational waves do not carry angular momentum. 

Notice that, with the assumptions we have made, we do not even know 
if w"' is a continuous vector field. We start with the final part of this sub- 
section, namely, to show that under a slightly stronger assumption on the 
differentiability of j"' on M, and on the metric qab at i, the tensor p"-^ given 
by (|TD satisfies Eq. (H). We have assumed that qab G W^'P^M). We now 



impose on the metric an extra condition given in Eq. (^. Then, we have the 
following result: 

Theorem 6. (Regularity on M) Assume that hypothesis in Theorem \^ 
hold. Assume that the metric satisfies If f G L'i{M), s"^ e W^''^{M), 
where q > 3, then, w°' G C^'"{M) and the tensor p"-'' satisfies PabV°^ll^ G 
L2(M). 

That w°' G C^'"(M) is deduced from standard elliptic regularity theorems. 
The second part is more difficult. The problem is that (Cqw)"' is not contin- 
uous at i, and so conditions that involve products of tensors are difficult to 
prove. Since the origin of the discontinuity in (Cqw)"' is the singular behavior 
of p"''', which we know explicitly, we proceed as follows. We split w"" into a 
regular part at i (called u°' in the proof) plus some divergent terms. These 
divergent terms are explicitly computed in terms of p"^ by an integration 
procedure based on Meyers' result p2[. Finally we show that this tu" satisfies 



a linear elliptic equation with a source in L'^{B^) with g > 3. Therefore it is 
C^''^ at i. Once this splitting near i on w"' is established, condition ( ^3|) is 
proved by explicit computation. 
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Proof. The source G L'^{M), while p"'', being smooth on M, belongs to 
L]^^j,(M). Therefore, standard elliptic regularity theorems given in [jl6[ imply 
e Wioc (M) c C^'°(M). 

We now begin the proof of the second part of the theorem. We work 
always with the rescaled metric qab and its corresponding covariant derivative 
Da- We explicitly compute the divergent terms of w"" at i. These terms are 
appropriate Meyers' potentials of the divergent terms present on Dijf^^y (See 
Lemma ^ below.) 

Let i?e C M be an open ball centered at i of geodesic radius e > 0. Let 
us choose e small enough such that the metric has the form (|31| ) in Riemann 
normal coordinates at i, and the cut function x is identically equal to 1 in 
B^. An explicit computation shows 

3 
k=2 

with = 0(1) and smooth on \ {z}, and continuous at i. The p^^^^ are 
functions of fia- We adopt the convention that a small circle over a quantity 
means that this quantity depends smoothly on fia, and does not depend on 
f. 

Let V^*^-) denote the Meyers potentials of p|j.^/f*^'^~^\ for each = 2,3, 

that is, vector fields V^*^-) = (^X]z^=o[^'^('^)]'''^(fco) /^^^"^^ defined on B^, with 
■y^^j^ appropriate functions of ha that can be explicitly computed in terms of 
p^Q^Y and satisfying 

So, here is our decomposition of the vector field w*, on B^, 

3 
k=2 



The rest of the proof shows that is indeed different iable at i. 
Thus a;* satisfies, 

3 

k=2 
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where (L^V(fc))* := (L^V(fc))^ - (LsVi^k))'- One can check that (L^V(fc))* = 
Q^^-{k-3)y Therefore the terms (L^V(fe))* with k = 2,3, belong to L^(M) 
with g > 3. Standard elhptic regularity implies that c<j' G W'^''^{M) C 
Ci'"(M). 

We have proved that the solution G C^'"(M) has the following expres- 
sion in B^, 

fe=2 \/=0 / 

Notice that the conformal factor uq is smooth on M, then an explicit com- 
putation implies that p"* given by Eq. (^1]) satisfies Eq. (pSf). □ 

We present here the generalization of Meyers' result, used in the proof of 
Theorem |^. 

Lemma 1. (Meyers' potential for L^) Consider the manifold {M.^,6ab), 
with 6ab the flat metric, da the metric connection, and let (LsV)"" = dbd^V°' + 
d"'dbV^/3. Consider the equation 

e 

{LsVr = r'^-'J2Mr)rpU^) (46) 

1=0 

where i > is a fix integer, r the geodesic distance from an arbitrary point 
p G M^, and p'^^^{n) is a C^'"(]R^) function of ria = daV, with K >0. 

Then, there exists C^~^'^'"{M.^) functions V^^{n), with I = 0, ■■■,£ + 2, 
such that 

e+2 

V'^ = r''Y.Mr)f V(?)H (47) 

1=0 

is a solution of ^4^). 

Proof. We look for solutions of Eq. (|46|) of the form 

V = v''- -d^X 
4 
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with 



n) 



1=0 



Lemma 4 in implies that there exist C^^'^''^{E?) functions wf^x, with 



/ = 0, 



+ 1, such that 



n] 



satisfies the first equation above. Then, one can exphcitly compute daV"", 
and again Lemma 4 in ||2^ imphes that there exist C^'^^'°(M^) functions 
A(/), with / = 0, ■ ■ ■ , + 2, such that 



1+2 



A = 5^ [ln(r)r A(ol 



n] 



is a solution of the second equation above. Therefore, an explicit computation 
gives (13) • □ 



3.3 Local regularity 

Consider a solution 6*, p"^ of Eqs. (|lO|) -(p!^). Assume now that the free data 
s'^^, jl and are f2-piecewise and tangentially smooth. In the first part of this 
subsection we then prove that the fields 9 and p"* are f2-piecewise smooth. 
(Theorem j^.) This proof is based on standard elliptic regularity theorems. 
In the second part of this subsection we prove that these fields are also VL- 
tangentially smooth. This result is split into two parts, first for linear elliptic 
systems, (Lemma 0), and then for Eqs. (p!0|)-([T3|). (Theorem ^j.) 

Theorem 7. (f2-piecewise smooth) Let both, qab and s"-'', in C°°{M). Let 
9 and p"^ he solutions of (|l^)-(|l^) given by Theorem ^ and Theorem |^. // 
the source functions fi and are Q-piecewise smooth then so are 9 and p"'^ . 

Proof. By the assumption qah € C°°(M) we have that the two elliptic opera- 
tors Lf^ and Lg have smooth coefficients in M. Applying the standard interior 
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elliptic regularity to the domains Q and M \ f2 we obtain that if j"" is Q- 
piecewise smooth then w"" is also f2-piecewise smooth. Because p"'^ G C°°(M), 
and by assumption s"'' G C°°{M), then p"'' is f2-piecewise smooth. 

In the case of 6 we note first that by the elliptic regularity 7 is smooth in 



M. Consider now equation (|22|) for Denote by f{x,-&) the right hand side 
of this equation. By the assumption on fi and the previous argument regard- 
ing w"", we have that the function f{x,i)) satisfies the following property: if 
belongs to C^'"(fi), [or to C*'"(M\f2)] then the composition f{x,-^{x)) de- 
fines a function that belongs to C^'"(f2) [or to C*'"(M\f2), respectively]. By 
Theorem ^ we know that the solution -i? G C^'"(M). (The argument works 
also with '(9 G C°(M).) Then we make an iteration, applying the elliptic 
regularity for the domains Vl and M \ in each step, to obtain that is 
fi-piecewise smooth. Therefore, so is 6. □ 

Let Q G Q', and V"' a smooth vector field on Q'. Let u any ten- 
sor field on M. Denote V^°\u) := u, V^^^{u) := VDaU, and V^''\u) : = 
VDaiV^^-^Xu)], for A; > 1. In order to prove tangential regularity we prove 
first the following lemma. 

Lemma 2. Let L he a linear elliptic operator of second order on some open, 
hounded set Q' (Z M with smooth coefficients. Let V"" he a smooth vector field 
on Q's, with Q C Q' . Let u G W'^'^{VL'), with p > 1, he a tensor field on Q! 
solution of the elliptic equation L{u) = f . Let k >0, an integer. 

(i) IfV^''\f) G LP(fi'), then V'^^\u) G W^^^i^'). 

(11) IfV^'^^f) G C"(fi'), then V^''\u) G C^'^i^l'). 

Proof. The proof is by induction on k. Consider the part (i) of the Lemma. 
The case = is the standard interior elliptic regularity. See for second 
order elliptic equations and [0, 0, ^ for systems. Assume now that (i) is 



true for k — 1. Consider now the following identity 

V^'\L{^)) = E (^)^'nV^,^](^^'-'H^)), (48) 

where (^) = k\/[l\{k — /)!] and we have introduced the notation 
^'\V,L]{u) ■.= L{u) 

(i)[\/,L](n) := [V,L\{u) = V{L{u)) - L{V{u)) 
('+i)[\/,L](w) := [\/,«[\/,L]](n). 
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Notice that, for all / > 0, the operator ^'-'[VjL] is a second order operator 
with smooth coefficients on M. Assume now that V^^\f) G LP{n'), and 
V^^\u) e W^'P{n'), for all </< A; - 1. If we write the identity (|8D as 

= V^'\f) - J2 (f)^'nV,L]iV^''-^\u)) (49) 
1=1 ^ ^ 

then all the terms in the right hand side belongs to Lp{Q'). Then the elliptic 
regularity theorems imply that V^'^^u) G ly^'^(n'), p > 1. The case (ii) is 
similar. □ 



Theorem 8. (i7-tangentially smooth) Assume the hypothesis on Theo- 
rem [2- // /i and j"- are Q-tangentially smooth then so are the fields 9 and p""^ 
which solve Eqs. (|TU|)-(|T§). 

Proof. Fix any open set in M such that Vt C Vt' . Let V to be the tangent 
vector field Vqq^ defined in Sec. Y,.l\ Since Eq. (Q) is linear, Lemma ^implies 
that w° is f2-tangentially smooth. 

Eq. (|2^ ) is semi-linear. However, there exists a solution d G W'^''^{VL') 
for g > 3. Therefore, d G C^'"(n'). This is the subtle step. Because d G 
C^'"(fi'), it implies that V{f{x,^{x))) G L''{yL'), for g > 3. The reason is 
that when we compute V{f), appear terms of the form "function in LP(f2')" 
times "V"('i?)." If d was only continuous, then, these terms would be not, in 
general, in LPiVt'). Then Lemma | implies that 1/(^9) G W'^'^Vt') C C^'^^iyt'). 
Thus V'^'^\'d) G C^iVL'). Then V'^'^\f{x,'d{x))) G Li{yL') and we obtain 
V^'^\'d) G C^'"(f2'). Iterating this argument, the conclusion follows. □ 



4 Further requirements 
4.1 Energy condition 

In order to understand the origin of this discussion on energy conditions it 
is useful to compare the usual procedure to find solutions of the constraint 
equations with matter sources. In that procedure, one rescales both, the 
energy density fi as well the momentum current density The rescaled 
j'^ = O'^^j"" is fixed from the requirement that the momentum constraint be 
independent of 9, while the rescaled jl = 0~^ii is chosen such that j / jl = j/ fj,, 
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where we have introduced j := a/ qabj"'j^, and j := a/ Qabj'^j'^- Therefore, if 
the energy condition is satisfied by the rescaled fields j"" and /i, with respect 
to the unphysical metric, then the physical fields j° and fl do satisfy the 
energy condition. That is the usual procedure. Here we can not rescale the 
energy density, because we need the extra condition that the physical energy 
density /2 be constant at the border of its support. 

We show here that the lower bound for jl given in (iii) in Theorem |^ is 
sufficient to guarantee that the physical matter fields satisfy the dominant 
energy condition. The idea is this: Initial data with = and fi satisfy 
trivially the energy condition. Therefore the same applies for arbitrary small 
j. Condition (iii) is just a rough bound on the smallness on j that also guar- 
antees the energy condition. Let 7_ = min-^i'j), with 7 the Green function 
solution of Eqs. (120)- (BT]). Then we have the following: 



Lemma 3. Let M , M , qab, /i, and p"-^ he as in Theorem 0. Let 9 he the 
corresponding solution of Eqs. (|71|^, ([7^. Let jl and have support in Vl 



and in C^{VL). If j < po(7-)^ then the fields Qab = O^Qab, = ^'^j"', and fi, 
satisfy i < fl. 

Proof. Since d is positive (see Theorem ^ > 7 on M. Then we have 

~3 = O-'j < <Po<fi. 

□ 



4.2 Inversion of Eqs. (|5|)-(§) 

We show here that Eqs. (||)-(^ are invertible, that is given the functions 
/i, then there exists unique functions p, satisfying these equations. In 
other words, the fluid 4-momentum density as seen by an arbitrary observer 
determines the fluid co-moving 4-momentum density. It turns out that the 
proof is not obvious and we did not find it in the literature. 

The main difficulty is that the map defined by Eqs. (i|)-(i|) is non-linear. 
Furthermore, it contains an unknown function, the state function, subject to 
minimally restrictive properties. Since v"" and j° are parallel, these equations 
reduce to 

i = (51) 
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with j = ^/ jaj"" and v = y/vaV"-. We define the map $ between subsets of 



^2 as 



Eqs. (|50|)-(|5TD can be rewritten as (/i, j) = $(p, v). Given a positive constant 
Po we define the following two subsets of 

D : = {(p, ^) G : Po < P, < S < 1} (53) 
/: = {(pJ)eM2:po(j)<p, 0<j}, (54) 



where po(j) '■= Po/2 + y Po/4 + j^- Our result is: 

Theorem 9. Let p{p) be a state function such that: (i) p{p) > for 
P > Po > 0; (ii) p(po) = 0; (in) < dp/dp < 1. Then, the map ^ : D ^ I 
is a difjeomorphism. 

Proof. First we prove that $ is bijective. 

Surjectivity: The essential tool is Brouwer's fixed-point theorem. (See for 
example [^.) Eqs. (^)-(^) are equivalent to 



p = fl-jv (55) 

(56) 



j p + pv"^ 



p p + p 

Fix a point (p, j) G /. Consider the map F given by 

' j X + p{x)y'^^ 



F{x,y) :-- 



(p-jy), 



p X + p{x) 



Introduce the compact convex set C := [po,p] x [0, 1] C M^. 

We claim that F -.C ^C. We write F{x, y) = (Fi(x, y), F2{x, y)). Then, 
by definition of /, j/p < 1 and so, for all (x, y) G D, we have < F2{x, y) < 1. 
We now show that, for all {x,y) E D, po < Fi{x,y) < p. The assumption 
< y < 1 implies p — j < p — jy < p. But po < Po{j) — j < p — j- Therefore 
po< p-jy < P- 

The map F is also continuous. Therefore, by Brouwer's fixed-point the- 
orem, there exists a fix point F{x,y) = {x,y). 
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Notice that j < flo{j) ^ implies that there exists no fix point of the 
form (x, 1). Therefore, we conclude that, given a point (/i, j) G /, there exists 
a point {p,v) G D which solves Eqs. (pOD-(pl|). 

Injectivity: Consider Eqs. (0)-(0), written as 

fj' = P + jv 
{p + p)v 



J 



1 — V 



1 )2 



Assume that there exist two points (pi, Vi) and (p2, "^2) which are solutions of 
these above equations for the same value of (p, j). If ■Ui = then the second 
equation above implies j = 0, and so V2 = 0, which in turn implies pi = p2- 
If vi = V2 then the first equation below implies pi = p2- 
Assume now that 7^ 0, ?)2 7^ 0, and V2. Then 

(P2 - Pi) +i{,V2 - = 



(P2 -Pi)(l -T^)= 3 



V2 Vi 



where z/^ = {dp/dp)\pi, with p' G [pi,P2] and we have used the mean value 
theorem for p{p) . Then the above equations and the assumptions on vi and 
V2 imply 

u'^ViV2 = 1- 

But by assumption z/^ < 1, so that we have a contradiction. Therefore, 
injectivity follows. 

It remains to prove that $ and are differentiable. By direct compu- 
tation and the assumptions on p{p) one can check that the derivative map 
of $ is invertible at each point of D. Then, by the inverse function theorem, 
is also differentiable. □ 

Notice that the proof fails if po = because the derivative of $ is not 
invertible at this point. This will be the case for an equation of state of 
the form p = ap"' , where a and 7 are constants. On the other hand, in this 
work we are interested in equations of state of liquid-type, i.e. such that the 
pressure vanish for a positive value of the density at the border of the fluid. 
For example p = a[(p/po)^ — 1] • For suitable constants a and po this equation 
describes water. (See |p!0| .) 
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5 Discussion 



The principal interest in the initial data given here is to use them to set 
up an initial value formulation. This formulation should be able to describe 
isolated, nearly static fluid bodies. That is why we have concentrated on 
finding liquid-type data and, inside this class, the smoothest possible data, 
i.e. the simplest to evolve. We have shown here that these data are not 
simple to obtain. 

The discontinuity of the fluid energy density at the boundary of its sup- 
port and extra constraints at that boundary [see Eq. (|^)] were the main 
difficulties. One main idea was to not rescale the fluid energy density and 
so, being free data, trivially solve the extra constraint Eq. (J^). (While also 
requires that the fluid 3- velocity vanish at the body boundary.) This un- 
conventional rescaling of the fluid flelds introduces difficulties in the task of 
finding solutions to the Hamiltonian constraint. These difficulties were solved 
in Theorem |^. Smoothest liquid-type data are almost-smooth, i.e. smooth 
except in the normal direction to the body boundary. The main step in 
establishing this result is Lemma ^ The rest is standard elliptic regularity. 

We have shown that at the body boundary, the first fundamental form is 
only in Wi^^{M), q > 3. This differentiability is below the threshold required 
by the known existence theorems on symmetric hyperbolic equations to prove 
existence of solutions associated with such a data. Theorems in Sec. 5.1 



m 



30| require initial data in with s > 5/2, where fl' C 



p3 



IS 



open and bounded. (See for a related improvement of this result and 
also [|18| for a discussion on the possible future development of the subject) 
Imbedding Theorems imply that W'^i^Q') C W^'P{Q'), with p > 3, but not 
in the other way around. Therefore, data in W'^'P{Q') is not enough for the 
known theorems to guarantee existence of solutions. 

We guess two possible ways to set up an initial value formulation for 
these liquid-type bodies. The first one is to study in detail the Einstein- 
Euler system given in |T^, with the hope that particular features of this 
system allows enough decrease in the differentiability threshold on the initial 
data to include the data given here. A second way is to set up two initial 
boundary value formulations, one for the interior of the body and one for 
the exterior, and then match both, in an appropriate way, at the boundary 



of the body (see and [^). It is far from being clear if either of these 
alternatives works. 
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A Appendix 

We discuss here the bound on the physical energy density given by Eq. (|^) in 
Theorem ^. In the first subsection, we show an inequality that is true for all 
maximal, asymptotically flat initial data with matter sources. This inequality 
is similar to Eq. (^) in the sense that relates the same quantities, but only 
the L^{Q) norm of the energy density appears. In the second subsection we 
show that ( plD is in fact a restriction, that is, there exist solutions of the 
constraint equations which do not satisfy it. Nevertheless we give arguments 
to show that physical systems like neutron stars do satisfy the bound (0). 

A.l An inequality 

Consider the following result. 



Lemma 4. Let M , M , and qab be as in Sec. HJ. Let p"^ and 9 he any 
solution of (|7^-([7^^ with pa"" = 0. Fix a point p e M, and denote by Br an 
open ball centered at p, of geodesic radius r. Then, for a sufficient small r, 
we have 

46 27rr 

MlHb.) < (57) 
where 7_ := infoB^'J, where 7 is defined in Eq. 
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Proof. Consider any solution, 6 of Eq. ([TT]) on M. Then, 



r, < ^^^^^ r^si 



Introduce as in Subsection |3.1| , that is, ^ = 7 + t?. We parametrize all 
possible solutions, instead by ^, by a function a := Lg{{}). Denoting by 
L^^(cr)(x) = — (l/47r) J^^a{y)-y{x — y)dV{y), then, inequality above trans- 
lates into 



K (J 



[7 + ^,-n^)P' 

The Green function has the form, 7(x — y) = l/\x — y\ + g, where g > 
on Br- (See pO|, ^.) The inequality 7(x — y) > l/(2r), that is true for all 
x,y & Br, implies that L^^{(t) > ||cr||ii(B^)/(87rr), and this, in turn, implies 

4"^""^^^^^ - [7- + lklUMB.)/(8vrr)]5- 

The last step of the proof is to maximize the right hand side of inequality 
above with respect to all possible functions a. The maximum value is taken 
for ||cr||Li(B,,) = 27rr7_, and inequality above gives Eq. (|57D. □ 

A. 2 Static spherical body 

In this subsection we explicitly construct an initial data set for a static, 
spherically symmetric, liquid-type body. We match, in appropriate coordi- 
nates and in a way, a 3-sphere endowed with its standard metric, with a 
3-dimensional Schwarzschild slice. The reason for re-doing this known con- 
struction (see [^) is twofold. First, this example, for suitable choices of 
the parameters, violates the bound (|2^). Second, we want to answer the 
following question: What kind of physical systems satisfy the bound (U)? 



We show that the answer turns out to be (at least for this example): stars 
with radius R > l.OSRg, where Rg = 2m is the Schwarzschild radius and m 
is the total mass. Note that this bound is bellow to i? > which is the 
necessary condition for hydrostatic equilibrium in General Relativity (see for 
example then this bound is expected to be satisfied for every star near 
equilibrium. 
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Let M = S^, the conformal metric be the standard metric, unit radius, 
of S^, the point i be the North Pole of S^, and the domain f2 be a ball centered 
at the South Pole of S^. 

Let 6ab be the flat metric, and r be the corresponding spherical radius. 
Consider the following initial data set 

qab = e%t, P'"' = 0. (59) 

The conformal factor 9 is given by 

where the positive constants a, tq, ri, and m, satisfy the following relations 

2 2rl (ro + m/2)6 



m ' 2mrQ 



(61) 



Using (|6l[) one check that 6', and hence ga^, is a function in M.^. There 
are two free parameters, for example we can take m and ro, m is the total 
mass of the data. The metric qab for r > ro is the Schwarzschild metric in 
isotropic coordinates and for r < ro is the standard metric on of radius 
a^/^. The Ricci scalar of the metric Qab is given by 

«=r. 'IT"' m 

if r < ro, 
and the physical energy density is 

if r > ro, 

f^= \ 6mrg .r ^ (63) 

K(ro+m/2)6 r <- ro. 

We see that the energy density fi has support in a closed ball of radius ro. 
In order to make contact with the assumptions in theorem § we write this 
initial data as follows. Using the well known relation 

7Vab = Sab, 7=P^ , (64) 



34 



we obtain 



For convenience, we have chosen a different normahzation for the Green 
function 7 than ( ]20|) and (|21|), in order to fix to be exactly the unit radius 
standard metric of S^. This difference in the normahzation wiU, of course, 
play no role in what follows. 

We want to prove that for some choices of the free parameters tq and 
m, this initial data violate the bound (PH). In order to do that we calculate 
explicitly the right and left hand side of (0). Since ji is constant we have 

m\mn) = mo\.(Sl)Y'\ (65) 

where Volgo(f2) denotes the volume with respect to the metric From Eq. 
(|6^ ) we have 



Since Cp = in this example, using ( |65D and ([66| ) we obtain that inequality 
([2^ ) is equivalent to 



[\o\o{n)Y"<P{l + ^)\ (67) 

where 13 = 4^/(3- 5^k) ^ 0.77 (the constant k, which depends only on 5*^ and 
can be calculated explicitly for this case k = \pX). Note that Volgo(f2) 
depends only on the dimensionless parameter mo/(2ro). For mo/(2ro) = 
we have that Volqo(f2) = Volgo(S'^) = 2-7?^ > then there exist values of 
mo/(2ro) such that the bound ( pi] ) is not satisfied. We use that Volqo(i7) < 
Volgo(5'^) for arbitrary mo/(2ro), to obtain a sufficient condition in order to 
satisfy Eq. (|67|) 

W(2ro) > 1-75. (68) 

Since the exterior metric is the Schwarzschild metric with mass m we can 
write this condition in terms of the physical radial area coordinate R to 
obtain 

i?> 2.16m. (69) 
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